AD-786  647 

EXPERIMENTAL  PULSED  LASER,  REMOTE 
CROSSWIND  MEASUREMENT  SYSTEM  - 
FEASIBILITY  STUDY  AND  DESIGN 

J.  Fred  Holmes,  et  al 

Oregon  Graduate  Center 


Prepared  for: 

Army  Electronics  Command 


July  1974 


DISTRIBUTED  BY: 


National  Technical  information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Roa'i,  Springfield  Va.  22151 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (1W)en  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

3.  RECIPIENT'S  CATALOG  NUMBER 

ht>  7<PS  £47 

«.  TITLE  (and  Subtitle) 

PULSED  LASER,  REMOTE  CROSSWIHD  MEASUREMENT  SYSTEM 
Feasibility  Study  and  Design 

S.  TYPE  OF  REPORT  A  PERIOD  COVERED 

FINAL;  7  December  1973  - 
30  June  1974 

8.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUI  HORfaJ 

J.  Fred  Holmes  and  J.  Richard  Kerr 

».  CONTRACT  OR  GRANT  NUMBERf#; 

DAAD07-74-C-0094 

9.  performing  organization  name  and  address 

The  Oregon  Graduate  Center 

Department  of  Applied  Physics  &  Electronic  Scieno 

t0.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  6  WORK  UNIT  HUMBERT 

1 1.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

U.S.  Army  Electronics  Command,  Atmospheric 

Sciences  Laboratory,  White  Sands  Missile  Range, 
New  Mexico 

12.  REPORT  DATE 

July  1974 

IS.  NUMBER  OF  PAGES 

42 

M.  MONITORING  AGENCY  name  &  AOORESSfi/  dllletent  from  Controlling  Ot/lca) 

i _ 

15.  SECURITY  CLASS,  (ot  thia  report) 

Unclassified 

15*.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 

16.  DISTRIBUTION  S.  XTEMENT  (ol  thia  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (o i  the  abstract  entered  In  Block  20,  It  different  from  Report) 

ie  supplementary  notes 

Detail*  of  illustration!  to* 
this  document  may  be  better 
studied  on  microfiche 

19  KEY  WORDS  (Continue  on  reverse  aid •  If  necessary  end  identify  by  block  number) 

Remote  Crosswind  Measurement,  Laser  Anemometer,  Atmospheric  Turbulence 

.  '  ‘  '  i  *  >  ’  f  i  , 

20.  ABSTRACT  f Continue  on  reverse  aide  If  necessary  end  Identify  by  block  number) 

The  feasibility  determination  and  design  for  an  experimental,  pulsed  laser, 
diffuse  targe1;,  remote  crosswind  measurement  system  is  described.  The 
systesn  consists  cf  a  laser  transmitter,  a  diffuse  target  and  a  receiver. 

a  unique  scheme  utilizing  a  double-pulsed,  Q-switched  ruby  laser  is  used  to 
measure  the  slope  of  the  time  delayed  autocovariance  function  at  zero  time 
delay  which  under  appropriate  conditions  is  proportional  to  the  weighted 

DD  I  IAN  7?  >473  ED.TION  OF  1  NOV  65  IS  OBSOLETE  /  „  ,  ,  , 

J  '  Unclassified 


r.r  '  JRiTY  CLASSIFICATION  OF  THIS  PAGE  flWicn  Data  Enteral; 


_ Unclassified _ 

SECURITY  CLASSIFICATION  OP  This  PAGCfSKSwi  Data  EntwmQ 


/ 


20.  Abstract  (continued) 

average  crossvind  along  the  laser  propagation  path.  The  receiver  consists  of 
two  silicon  vldicons,  a  spinning  disk  for  directing  alternate  pulses  to  the 
vid icons  and  appropriate  scanning  circuitry  for  generating  a  10  x  10  array 
of  apertures  on  each  vidicon.  Due  to  the  pulsed  nature  of  the  system, 
spatial  (ensemble)  averages  in  lieu  of  time  averages  are  utilized.  The 
system  can  measure  the  autocovariance  function  and  the  slope  of  the  time 
delayed  autocovariance  function  at  zero  time  delay  for  nine  aperture 
spacings  in  two  orthogonal  directions  and  has  the  potential  for 
two-dimensional,  path  resolved  crossvind  measurement . 


j 

i 

i 


•  > 


...  __lto.cla.ss  if  led _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGEITrtion  Dalit  Enlnrd) 


Experimental  Pulsed  Laser,  Remote  Crosswlnd 
Measurement  System  —  Feasibility  Study  and  Design 

Dr.  J.  "red  Holmes 
Dr.  J.  Richard  Kerr 


Contractor:  Oregon  Graduate  Center  for  Study  and  Research 

Sponsor:  U.  S.  Army  Electronics  Command,  Atmospheric  Sciences 

Laboratory,  White  Sands  Missile  Range,  New  Mexico  38002 

Contract  Number:  DaA  D07-74-C-0094 

Effective  Date  of  Contract:  7  December  197  ; 

Contract  Expiration  Date:  30  June  1974 

Amount  of  Contract.:  $65,000 

Principal  Investigator'  Dr.  J.  Richard  Kerr 

Co-Investigator:  Dr.  J.  Fred  Holmes 

Project  Monitor:  Mr.  Thomas  H.  Pries 


Approved  for  Public  Release; 
Distribution  Unlimited 


% 

i* 


v 


CONTENTS 

Page 


INTRODUCTION .  1 

BACKGROUND .  3 

EXPERIMENTAL  PROTOTYPE  DESIGN  .  9 

PROCESSING . 23 

EXPERIMENTAL  PROGRAM . 28 

REFERENCES . 36 

APPENDIX  A . 38 

FIGURE'?  1-12 

Figure  1  -  Signal  to  Noise  Ratio  Versus  Gate  Time  .....  13 

Figure  2  -  Receiver  Diagram  .  16 

Figure  3  -  Transmitter  Diagram . 18 

Figure  t  -  Receiver . 20 

Figure  5  -  Microprocessor . 21 

Figure  6  -  Transmitter  Installation  .  ...  22 

Figure  7  -  Field  Site  Facility . j0 

Figure  8  -  Calculated  Covariance  Curve . 31 

Figure  9  -  Calculated  Wind  Weighting  Function  .......  32 

Figure  10  -  Scintillation  Pattern,  500  Meters  .  33 

Figure  11  -  Scintillation  Pattern,  500  Meters  .  34 

Figure  12  -  Scintillation  Pattern,  1.6  Kilometers  .....  35 


iV 


TNTRODUCTION 


The  feasibility  of  remotely  determining  the  transverse  wind 

12 

velocity  with  a  point-to-point,  cw  laser  system  is  well  established.  ’  * 
4  5 

In  recent  work,  ’  investigators  have  considered  the  extension  of  the 
technique  to  the  use  of  a  cw  laser  as  well  as  ambient  light,  with  a 
passive,  diffuse  target.  In  either  case,  the  technique  basically  con¬ 
sists  of  relating  the  scintillation  pattern  on  the  receiver,  to  the 
frozen-in,  transverse  motion  of  the  turbulence  structure  along  the 
propagatioi,  path.  Tho  most  successful  cw  method  that  has  been  used 
involves  tne  measuremen  of  the  slope  of  the  log  amplitude  covariance 
function  at  zero  time  de.'ay.^ 

The  field  application  for  the  wind  measuring  technique  requires 
nearly  instantaneous  measurements.  Although  short  time  averages  (e.g. 
1-10  seconds)  have  been  shown  to  produce  satisfactory  results  in  the  cw 
approach,  the  use  of  a  non-cooperative,  diffuse  target  suggests  the 
utilization  of  a  pulsed  laser  cource  for  useable  signal  levels  at 
reasonable  ranges. 

In  order  to  employ  a  pulsed  laser  and  still  make  a  temporal  measure 
tnent  as  required  for  wind  velocity,  a  double-pulse  technique  is  used 
whereby  two  Q-switched  laser  pulses  are  generated  at  a  time  interval 
on  the  order  of  a  millisecond.  This  interval  appears  to  be  appropriate 
for  the  employment  of  the  covariance  "slope  at  zero  time  lag"  technique"*" 
and  is  compatible  with  current  double-pulse  laser  technology. 

Since  the  pulse  technique  precludes  time  averaging  of  the  statis¬ 
tical  quantities  measured,  spatial  averaging  is  employed.  Therefore 
the  equivalent  of  an  anay  of  detecr.ors  had  to  be  utilized  in  order 
to  generate  a  sufficient  number  of  sample,;  of  the  instantaneous 
scintillation  pattern  to  use  spatial  averaging.  To  the  extent  that 
the  statistics  of  this  pattern  are  invariant  with  lateral  translation 
at  the  receiver  (a  safe  assumption  with  the  use  of  a  diffuse  target 
as  a  secondary  source),  such  spatial  averaging  wiL]  give  the  desired 
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result  and  the  pulsed  technique  is  fundamentally  similar  to  the  cw, 

time-averaging  approach.  However,  the  partial  coherence  of  the  laser 

and  the  characteristics  of  scintillations  involving  a  diffuse  source 

introduce  new  factors.  These  issues  are  pertinent  foi  both  the  cw 

and  pulsed  systems  and  have  been  given  some  consideration  by  other 
4  5 

investigators.  * 

The  primary  emphasis  during  Phase  I  cf  this  investigation  was  on 
feasibility  determination  and  design  of  an  experimental,  pulsed-laser , 
remote-crosswind  measurement  system  utilizing  the  "slope"  method. 
Accomplishments  during  the  report  period  include: 

Established  a  receiver  concept  for  real  time  measurement 
of  the  average  wind  utilizing  a  pulsed  source  and  diffuse  target; 

Established  that  use  of  a  diffuse  target  with  pulsed  laser 
illumination  for  wind  determination  is  feasible; 

Designed  and  built  an  experimental  transmitter  and  receiver 
system; 

Established  a  test  site  and  installed  and  tested  the  trans¬ 
mitter; 

Developed  a  technique  for  processing  the  data  to  determine 
the  wind;  and 

Established  a  conceptual  design  for  a  field  system  requiring 
realistic  laser  energy  levels. 
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BACKGROUND 


In  r.cn-saturated  turbulence,  the  time  delayed  covariance  function 
t„>r  the  intensity  is  given  by^’^*^’^ 


Cy  (P,t)  =  4tt2  0.033k2  f  C2  (z)  W2  (P,t,z)  dz 

~  c\ 
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The  parameters  used  in  equations  (1)  -  (4)  are  defined  by 
Kj^2  =  Receiver  Filter  function 
Kg2  =  Source  filter  function 
Dr  =  Receiver  aperture  diameter 
Dg  =  Incoherent  source  diameter 
L  =  Total  propagation  path  length 
z  =  Distance  from  source 
Cn2(z)  =  Turbulence  strength 
V(z)  =  Wind  velocity 
p  =  Receiver  aperture  separation 
t  =  Time  delay 
k  -  Wave  number 

A  Kolmogorov  spectrum  has  been  assumed. 

The  slope  with  respect  to  t  of  the  time  delayed  covariance  function 
can  be;  found  from  (1)  and  is  given  by 


M  (p)  -  3CX(p>t) 
3t  ~ 


T  -  O 


(5) 


My  (p)  ■  (0.033)  4Tr2k2  j  v(z)Cn2  (z)W2  (z,p)  dz 

J  o 
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where 
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Equations  (1)  and  (6)  hold  for  the  log-amplitude.  However,  it 

is  the  intensity  that  is  measured.  This  can  be  accounted  for  by 

6  8 

normalizing  (1)  and  (6)  by  the  variance  since  it  has  been  shown  ’ 
that  the  normalized  covariance  of  the  log-amplitude,  amplitude  and 
intensity  are  equal  for  the  non-saturated  casev  The  normalized 
expressions  are  then  given  by 


CTM  (p',t)  =  C_(p',x)  -  ;  C2(W)F!(W,P  ,t>  dw 
IN  1  Jo  n 
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and 


V2J> 


(12) 


and  where  the  following  variable  changes  have  been  made 
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Equations  (7)  and  (8)  represent  the  theoretical  basis  upon  which  remote 

probing  of  the  turbulence  and  crosswind  velocity  is  based.  Equation  (7) 

could  be  used  to  determine  the  average  wind  velocity  by  finding  the  time 
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delay  x  at  which  it  has  a  peak  value.  ’  However,  if  the  turbulence 
and/or  wind  is  non-uniform,  multiple  peaks  occur.  Also,  the  scintillation 
pattern  tends  to  evolve  with  time.  Consequently  this  method  is  not  very 
useful . 

Equation  (7)  is  linear  in  C  2 (a)  and  non-linear  in  V(z).  It  could 

n  9 

be  used  to  determine  the  wind  velocity  ilong  the  path,  but  the  non-linear 
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relationship  to  V(z)  makes  it  difficult.  At  zero  time  delay,  equation 

(7)  is  independent  of  the  velocity  and  depends  linearly  on  C  2(z). 

n 

The  "slope"  function,  Equation  (8),  is  linearly  dependent  on  the 
product  of  wind  velocity  V(z)  and  the  turbulence  Cn2(z).  Equations 
(7),  with  zero  time  delay,  and  (8)  can  at  least  in  principle  be 
solved  for  the  turbulence  and  wind  profiles.^  ^ 

If  the  turbulence  is  assumed  to  be  uniform,  then  C  2(z)  can  be 
taken  out  of  the  integrals  in  (8)  and  it  becomes 


(p  )  =  r1 


i 

V(W)F2(W,p'  )  dw 


1 

F3(W)  dw 


(13) 


The  "slope"  function  in  this  case  is  proportional  to  a  path  weighted 
average  of  the  velocity  distribution  along  the  propagation  path  and 
has  been  used  quite  successfully  for  cw  remote  wind  sensing.  ’ If 
the  wind  velocity  is  also  constant  along  the  path,  then  (13)  becomes 


M  (o'  )  =  6-1  V  F2(W,p'  )  dw  (14) 

N  _ 

..  1 

F3(W)  dw 
^  o 

The  slope  function  in  this  case  is  proportional  to  the  crosswind 
velocity. 

In  a  field  situation,  both  the  turbulence  and  the  wind  may 
be  non-uniform  along  the  propagation  path.  Consequently  a  field 
system  may  require  the  capability  for  path  resolved  wind  measure¬ 
ments.  There  are  several  possible  techniques  that  might  be  used 
to  accomplish  this.  In  principal  it  could  be  done  by  varying  the 
path  weighting  function  »^»16,  ^  the  weighting  function  takes 

on  the  form  of  a  "delta"  function,  then  the  turbulence  strength  and  wind 
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c 


velocity  would  be  determined  at  the  point  along  the  path  that  the 
weighting  function  was  non-zero. 

Another  approach  that  looks  very  promising  and  has  been  used 
9-11 

with  some  success  is  to  use  the  covariance  function  and  the 

slope  of  the  covariance  to  obtain  the  turbulence  strength  and  wind 
speed  as  a  function  of  position  along  the  propagation  path.  Each 
of  these  functions  is  represented  by  an  integral  involving  the 
wind  speed  and/or  the  turbulence  strength.  By  measuring  the  co- 
variance  and  "slope"  at  several  detector  spacings,  the  integral 
equations  can  be  used  to  resolve  C^(z)  and  V(z)  to  some  degree  along 
the  path.  An  important  issue  in  any  technique  is  the  required 
accuracy  on  the  measurements  versus  the  accuracy  of  the  path  resolved 
parameter. ^»11>1^»1®  As  the  path  resolution  increases,  the  required 
accuracy  of  the  measurements  increases  very  rapidly.  Consequently,  a 
practical  limit  appears  to  be  the  utilization  of  measurements  at  three 
or  four  detector  spacings. 
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EXPERIMENTAL  PROTOTYPE  DESIGN 


As  shown  in  the  previous  section,  the  covariance  function  and  the 
slope  of  the  time  delayed  covariance  function  at  zero  time  delay  contain 
information  about  the  strength  of  turbulence  and  the  crosswind  Velocity 
along  the  propagation  path.  Consequently,  an  instrument  that  can 
measure  these  two  functions  has  the  capability  for  remotely  determining 
the  turbulence  strength  and  crosswind.  This  is  the  basis  on  which  an 
experimental  prototype  was  designed. 

Work  by  other  investigators  has  been  restricted  to  cw  sources  where 

the  propagation  path  is  direct  to  the  receiving  system  (double-ended 

19  20 

system)  or  a  cooperative  target  such  as  retro-reflector  *  is  used 
(single-ended  system),  or  passive  (ambient)  illumination  is  used.  Under 
these  conditions,  an  acceptable  signal  to  noise  ratio  can  be  achieved 
at  reasonable  transmitter  power  levels.  However,  the  proposed  appli¬ 
cations  of  the  work  sponsored  under  this  project  requires  a  single- 
ended  system  operating  against  a  non-cooperative  target  which  can  be 
viewed  as  a  perfectly  diffuse  surface  with  reflectivity  of  approxi¬ 
mately  0.1.  Under  this  constraint,  a  cw  system  would  require  exces¬ 
sive  transmitter  power  levels.  However,  preliminary  nnalysis  indicates 
that  a  pulsed  system  operating  at  about  the  1  Joule  level  could  be  de¬ 
signed  with  adequate  signal  to  noise  'atio  (S/N). 

The  primary  noise  sources  in  the  system  are  the  Poisson  fluctu¬ 
ations  of  detector  current  and  thermal  noise  in  the  electronics.  For 
a  given  pulse  energy,  as  the  pulse  width  decreases  the  instantaneous 
power  increases.  Consequently,  the  detector  signal  current  increases 
and  eventually  its  Poisson  fluctuation  becomes  the  dominant  noise 
source.  If  the  receiver  is  gated  synchronously  with  the  transmitted 
pulses,  at  this  point  it  would  be  signal  shot  noise  limited.  There¬ 
fore,  in  a  direct  detection  optical  system,  the  signal  to  noise  ratio 

for  a  given  total  energy  increases  as  the  pulse  duration  is  made 
21 

shorter.  There  is  of  course  a  limit  on  how  short  the  transmitted 
pulses  or  ultimately  the  receiver  gate  width  can  be  made.  Tp  our 
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particular  case,  the  minimum  transmitted  pulse  width  is  around  10- "iO 

nanoseconds.  However,  it  will  not  be  necessary  to  go  to  this  limit 

in  gating  the  receiver  in  order  to  meet  the  system  requirements. 

For  a  photodiode  type  detector  working  into  an  ideal  (infinite 

bandwidth)  operational  amplifier  operating  in  the  current  mode,  the 
22 

S/N  is  given  approximately  by 


where 

en  ■  Amplifier  voltage  noise^/Hz 

e,  *  bias 
bn 

R^  *  Feedback  resistance 
r^  *  detector  shunt  resistance 
B  *  bandwidth 


R1  +  rd 

Rj  =  Amplifier  Input  resistance 

I  2  *  4KT  -  f  T  T  x 
n  —  +  2e  (Ig  +  Id  +  Ij) 

d 

I  =  detector  signal  current 
s 

1^  =  detector  dark  current 
Ij  =  amplifier  input  current 


voltage  nota&j/Az 
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If  I  2  is  dominated  by  2sl  and  R  2  I  2  dominates  the  other  de- 
n  J  s  f  n 

nominator  terms  in  (15),  then  the  system  operation  is  signal  shot  noise 
limited  and  (15)  becomes 


If  B  is  made  equal  to  1/t  where  t  is  the  pulse  width  (or  Gate  time), 

and  using  1 

u  en 

I  =  r - 

s  hvt 

where 

U  =  Received  energy 
R 

n  =  detector  efficiency 
equation  (16)  becomes 


S 


N 


V 

2hv 


(17) 


where 

U  =  U_  (Tarcet  Reflectance)  u  .  / Receiver  Diameter \  2 

R  T  4  l  1 

rt  \  Ran"e  / 

(18) 

For  a  target  reflectance  of  0.1,  a  receiver  diameter  of  2  cm  a 
transmitted  energy  of  1  Joule  and  a  range  of  3  Km,  signal  shot  noise 
limited  operation  yields  an  S/N  of  almost  60  db  which  should  be  ade¬ 
quate.  It  should  be  noted  that  (17)  depends  linearly  on  the  trans¬ 
mitted  energy.  Consequently,  if  the  transmitted  energy  is  decreased 
by  a  factor  of  ten  yielding  100  millijoules,  the  S/N  decreases  by 
only  10  db. 

Fie\.re  1  shows  the  effect  of  receiver  gate  time  on  the  signal 
to  noise  ratio  for  the  above  described  conditions  and  with  the  follow¬ 
ing  electronics  parameters  used  for  illustration 
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t  =  Gate  Time 


B  *  1/t 

=  106  ohms 

r ,  =  10s  ohms 
a 

Ri  =  rd 

I  =  U  /t 
s  r 

Id  =  5  x  10"8 

Ij  -y  Negligible 

e  =  30  x  10‘9  V-Hz"l/2 
n 

eL  -  30  x  10"9  V*Hz-1/2 
b 

n 


From  Figure  1  it  appears  that  signal  shot  noise  limited  operation 
can  be  achieved  with  a  gate  time  of  around  10“6  seconds. 

Several  types  of  detectors  were  considered  for  the  experimental 
prototype  including  silicon  vidicons,  arrays  of  photodiodes  and  charge 
coupled  devices.  The  presently  available  charge  coupled  devices  are 
several  orders  of  magnitude  inferior  from  a  noise  standpoint  than  the 
other  devices  and  were  rejected  on  that  basis.  The  silicon  vidicon 
(which  is  essentially  a  large  array  of  photodiodes  read  out  with  an 
electron  beam)  and  photodiodes  have  similar  noise  characteristics. 
However,  for  purposes  of  an  experimental  prototype,  both  have  advan¬ 
tages  and  disadvantages.  The  silicon  vidicon  is  limited  in  beam 
current  and  sweep  speed  which  in  turn  limits  the  minimum  gate  time 
and  prevents  signal  shot  noise  limited  operation.  However,  it  allows 
a  very  flexible  and  readily  changeable  design  due  to  the  need  for  only 
a  single  set  of  optics  (vice  an  array)  and  the  ability  to  electronically 
generate  and  control  the  receiving  array.  This  is  a  considerable  ad¬ 
vantage  and  led  us  to  choose  the  silicon  vidicoi.  for  the  first  experi¬ 
mental  prototype.  It  should  be  noted  that  at  this  point  it  appears  a 
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Figure  1.  Signal  to  Noise  Ratio  Versus  Gate  Time 


photodiode  array  will  be  necessary  for  a  field  system  due  to  the  need 
for  signal  shot  noise  limited  operation.  To  compensate  for  the  slower 
gate  time  (smaller  S/N)  with  the  silicon  vidicons,  a  special  target 
material  is  being  used.  It  is  a  high  gain  Scotchlite  paper  which 
represents  a  quasi-diffuse  scatterer.  When  illuminated  by  a  cw  co¬ 
herent  source,  the  Scotch1 ite  paper  exhibits  strong  speckle,  character¬ 
istic  of  a  diffuse  target.  Use  of  this  material  allows  the  effect  of 
a  diffuse  target  material  to  be.  assessed  at  a  much  lower  illumir.it ion 
level  than  would  normally  be  required  for  the  silicon  vidicon  system. 
The  Scotchlite  has  a  gain  relative  to  that  of  a  perfectly  diffuse 
surface  of  around  40  u.  The  signal  to  noise  ratio  for  the  silicon 
vidicon  system  is  estimated  to  be  in  excess  of  70  db  at  1.6  Km.  The 
maximum  background  illumination  that  is  expected  to  be  encountered 
at  the  tesc  site  is  0.03  Watt/cm2 /Sr/ym.  By  limiting  the  field  of 
view  of  the  receiver  to  3  milliradians  and  using  a  10“2  pm  bandpass 
filter  in  front  of  the  vidicons,  the  received  background  power  density 
is  of  the  order  of  3  x  10“ 10  Watts/cm2.  To  further  limit  the  back¬ 
ground,  a  mechanical  shutter  (spinning  disk)  is  used  to  mechanically 
gate  the  receiver  for  around  1  millisecond.  The  received  background 
energy  then  is  of  the  order  of  3  x  10“13  Joules/cm2  as  compared  to 
around  1Q“9  Joules/cm2  at  1.6  Km  for  the  received  signal. 

The  purpose  of  the  receiver  is  to  provide  intensity  measurements 
from  which  the  covariance  function  and  the  "slope"  of  the  covariance 
function  can  be  calculated.  The  latter  function  can  be  derived  by 
calculating  the  covariance  at  two  closely  spaced  points  in  time. 

This  is  accomplished  by  utilizing  two  transmitter  pulses  separated  in 
time  by  10“ 3  second. 

Since  the  quantity  being  measured  (received  intensity)  is  a  random 

variable,  a  single  measurement  would  not  yield  a  good  estimate  of  the 

expected  value.  With  a  cw  system,  a  time  average  can  be  used  to 

statistically  average  the  signal.  It  appears  that  good  results  are 

1  4 

being  obtained  with  temporal  averaging  times  of  1  to  10  seconds.  v 
Since  the  pulsed  system  samples  only  one  instance  in  time,  a  spatial 
average  must  be  used. 
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In  order  to  obtain  the  required  spatial  averaging,  thi.  equiva- 

5f> 

lent  of  an  array  of  detectors  is  required.  In  addition,  since  the 
two  pulses  are  only  a  millisecond  apart,  the  array  oust  be  read  out 
quickly  or  two  arrays  must  be  used  in  conjunction  with  a  means  for 
directing  the  return  from  the  second  pulse  to  the  alternate  array. 

We  are  using  two  silicon  vidicons  tc  realize  the  arrays.  The  face 
of  each  vidicon  is  divided  electronically  into  a  10  x  10  matrix  of 
equally  spaced  cells  Each  element  in  the  matrix  is  scanned  and  a 
digital  number  generated  that  represents  the  integrated  flux  over 
the  cell.  The  location  of  the  center  of  each  cell  (i.e.  the  spacing) 
is  fixed  in  the  present  arrangement  but  could  be  made  electronically 
variable.  However,  the  size  of  each  cell  is  controlled  electronically 
and  gives  us  the  capability  to  vary  the  effective  aperture  of  each 
cell.  A  spinning  disk  and  a  hole  and  mirror  arrangement  are  used  to 
direct  the  two  received  pulses  to  one  vidicon  or  the  other  (see 
Figure  2).  The  2  x  100  numbers  representing  the  two  scintillation 
patterns  are  temporarily  stored  in  a  memory  and  then  recorded  by 
punching  a  paper  tape  for  later  processing.  The  measurement  system 
is  extremely  flexible  which  will  aid  in  assessing  the  effects  of 
detector  spacing  and  aperture  size  on  system  performance.  It  will 
allow  a  variety  of  techniques  for  converting  the  information  in 
the  scintillation  patterns  into  wind  speed  to  be  evaluated  and  has 
the  capability  for  determining  two-dimensional  path  resolved  wind. 

The  receiver  is  illustrated  schematically  in  Figure  2. 

The  receiver  optical  system  consists  of  two  lenses  and  a  field 
of  view  limiting  aperture.  The  scintillation  pattern  at  the  entrance 
aperture  (first  lens)  is  demagnified  by  a  factor  of  approximately 
14.3  and  imaged  onto  a  1  cm  diameter  area  on  the  vidicon  face.  For 
a  1.6  Km  propagation  path,  this  yields  a  cell  spacing  and  cell 
aperture  of  0.43  Fesnel  zones  (v^AlT)  .  The  10  x  10  array  covers  an 
area  of  4.3  x  4.3  Fesnel  zones.  For  the  smallest  cell  spacing,  this 
allows  90  spatial  samples  to  be  averaged.  Of  course  all  the  samples 
are  not  completely  independent.  Data  from  individual  cells  can  be 
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combined  to  form  larger  cells  if  desired  and  the  second  lens  can  be 
changed  to  decrease  the  spacing.  The  primary  criterion  on  system 
performance  is  the  error  in  estimating  the  crosswind.  As  the  cell 
size  increases  more  signal  (average)  is  received,  but  the  scintil¬ 
lation  level  decreases.  Also,  for  a  given  cell  spacing,  as  the  cell 
size  increases,  the  samples  become  more  dependent  which  will  increase 
the  error.  Because  these  effects  are  competing,  there  is  probably 
an  optimum  cell  size  and  spacing  for  a  given  number  of  cells  and  the 
optimum  solution  may  also  be  a  function  of  range  and/or  turbulence 
level.  This  will  be  studied  during  phase  II. 

The  transmitter  is  shown  schematically  in  Figure  3.  The  source 
is  a  double-pulse,  Q-switched  ruby  laser.  The  requirement  for  the 
laser  to  produce  two  Q-switched  pulses  8' unrated  by  about  one  mini¬ 
s'  icond  was  met  by  exciting  the  lac?r  rod  with  a  long  pump  pulse  and 
Q-switching  twice  within  the  pump  pulse  interval.  A  double-pulse, 
Q-switched  ruby  laser  meeting  these  requirements  was  available 
commercially  (Apollo  Lasers  Inc.)  for  use  in  the  system.  The  energy 
per  pulse  is  approximately  0.5  Joule  with  a  repetition  rate  of  6  pulse 
pairs  per  minute. 

Ruby  appears  to  be  the  bast  choice  at  the  present  time.  However, 
in  the  future  longer  wavelengths  may  be  seriously  considered  because 
of  the  reduced  probability  of  saturation  of  scintillation  at  longer 
wavelengths.  A  small  He-Ne  laser  is  used  for  alignment  and  as  a 
tracer  beam  for  the  ruby  laser.  This  greatly  facilitates  steering 
the  beam  onto  the  target.  It  will  not  be  required  in  a  field  system. 

The  9.5  mm  output  beam  of  the  ruby  laser  is  beam  expanded  using 
a  negative  lens  and  a  320  mm  mirror.  This  reduces  the  3-5  milliradian 
divergence  of  the  raw  beam  to  around  10-4  radian  which  is  the  same 
order  of  magnitude  as  the  maximum  spread  due  to  turbulence.  A  nega¬ 
tive  l**ns  was  required  to  avoid  tightl>  focusing  the  raw  beam,  since 
it  has  enough  energy  to  break  down  the  air.  The  negative  lens  is  on 
on  ft  track  that  allows  the  beam  to  be  focused  to  a  spot  at  any 
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Figure  3.  Transmitter  Diagram 


distance  beyond  a  few  hundred  feet  from  the  transmitter.  The  spot  size 
in  the  absence  of  turbulence  should  be  around  160  mm  at  1.6  Km.  As  the 
turbulence  increases,  the  spot  si23  will  increase.  A  larger  spot  size 
decreases  the  degree  of  scintillation  and  changes  the  covariance 
function.  However,  the  larger  spot  size  is  caused  by  stronger  turbu¬ 
lence  which  increases  the  measured  effect.  Consequently  there  are 
competing  effects  which  tend  to  cancel  each  other.  The  change  in  the 
covariance  function  may  also  significantly  affect  the  instrument 
calibration. 

Another  factor  that  must  be  considered  is  saturation.  For  a 

reasonably  long,  low  path  at  visible  or  near-IR  wavelengths,  the 

scintillations  will  be  "saturated"  in  moderate  and  strong  turbulence.^ 

This  represents  the  multiple  scattering  realm,  and  affects  the  shape 

A  23 

of  the  covariance  curve  and  hence  path-weighting.  ’  For  a  system 

to  cope  with  this  phenomenon,  it  will  probably  be  necessary  that  the 

degree  of  sacuration  be  known  or  measured.  Our  measurement  system 

can  generate  an  estimate  of  the  covariance  function  at  9  spatial 

separations.  It  may  be  possible  that  the  degree  of  saturation  can 

then  be  inferred  from  the  covariance  function.  By  using  a  model 

24 

appropriate  for  the  saturated  condition  proper  corrections  can 

then  be  incorporated  into  the  data  processing.  In  this  regard,  the 

technique  for  determining  the  wind  from  the  data  may  have  to  be  self- 

25 

adaptive.  In  addition,  recent  work  by  NOAA  indicates  that  the 
effects  of  saturation  can  be  neutralized  by  appropriate  choice  of 
source  and  receiver  aperture  sizes.  More  study  is  needed  in  this 
area. 

Photographs  of  the  completed  equipment  are  shown  in  Figures  4-6. 
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Figure  5.  Microprocessor 
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Figure  6.  Transmitter  Installation. 


PROCESSING 


The  instrument  described  in  the  previous  section  records  200 
numbers  that  are  proportional  to  the  total  integrated  intensity 
received  by  each  cell  in  the  array.  From  these  data,  the  average 
intensity,  variance  and  covariance  function  for  each  pulse  and  the 
slope  of  the  covariance  function  at  zero  time  delay  are  calculated. 
The  covariance  of  a  random  variable  I(r,t)  can  be  expressed  as 

CI(r2,t2;  rlttj)  =  <[I(r2,t2)  -  I (r2 »t2)][I(r1,t1)  -  Krj  ,t  x ) ]>  (19) 

If  the  process  is  wide  sense  stationary  and  locally  homogeneous,  then 
the  covariance  function  only  depends  on  the  coordinate  differences 

p  =  |r2  -  rj | 
t  =  t2  -  tx 

and  (19)  becomes 

C].(p,t)  =  <[I(r2,t2)  -iHKrj.tj)  -  I]>  (20) 

The  data  is  reduced  on  the  basis  of  equation  (20). 

Since  we  are  dealing  with  .’■.wo,  two-dimensional  arrays  of  data,  it 
is  convenient  to  use  three  subscripts  to  identify  the  array  location 
and  pulse.  For  example  consider  I(J,K,N). 

J  =  array  row 

K  =  array  column 

N  =  pulse  number  (1  or  2) 

Using  this  notation  and  approximating  the  ensemble  average  by  a  finite 
spatial  average,  the  average  intensity,  variance  and  normalized  co- 
variance  function  are  then  given  by 
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where  S  is  the  normalized  cell  spacing,  n  takes  on  the  integer  values 
1  through  9. 

The  slope  of  the  covariance  function  at  zero  time  delay  can  be 
expressed  as 


M  (p')  = 
N 


9C  (f/.t) 
N 
at 


T  -  O 


=  limit  C  (ff.xo)  -  C  (p',-xo)  (24) 

x  -*  O  TJ 

0 - 

2xo 

With  2xo  chosen  to  be  10“ 3  seconds,  the  limit  restriction  can  be 
removed  without  causing  a  significant  error. *  Doing  this  and  rewriting 
the  covariance  functions  in  a  form  more  readily  identifiable  with  the 
arrays  and  pulses,  (24)  becomes 


M  (o')  =  c_  (X2,t2;X1,t1)  -  C  (X2,t1;X1,t2) 
N  _ 

2xo 


(25) 
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Now  using  (19)  this  becomes, 

MI  (p)  -  _±_  <[I(X2,t2)  -IJlKXLt!)  -I] 
N  2to 


-  [ I (X2 , t ! )  -  I][I(Xltt2)  -  I)>  /OjCN  =  l)a  (N  =  2) 

(26) 

Using  the  subscript  notation  and  approximating  the  ensemble  average  by 

a  spatial  average,  the  normalized  "slope"  function  becomes 

Mt  (ns*  ) 

N 

1q3  10-n  r  _ 

_ _ _  Z  Z  [I(J,K+n,2)  -  I(N  =  2)] 

10  (10-n)  a  (N  =  1)  a  (N  =  2)  J  =  1  K  =  1  L 
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[I(J,K,1)  -  I(N  =  1)]  -  [I(J,K+n, 1)  -  I(N  =  1)][I(J,K,2)  -  I(N  =  2)] 

(27  ) 


The  data  stored  in  each  array  is  not  the  intensity  but  a  number 
proportional  to  the  intensity.  In  addition,  the  proportionality  con¬ 
stant  will  be  different  for  each  array.  However,  because  of  the  way 
in  which  the  covariance  and  "slope"  functions  are  normalized,  the 
results  for  these  quantities  are  independent  of  the  proportionality 
constants  and  the  energy  in  each  laser  pulse  and  no  related  cali¬ 
bration  is  necessary.  However,  this  does  not  compensate  for  any 
sensitivity  variation  across  the  tube  face.  If  absolute  data  on 
the  average  intensity  and  variance  is  desired,  then  the  two  pro¬ 
portionality  constants  will  have  to  be  either  measured  or  calculated. 

To  obtain  the  average  wind,  the  measured  slope  is  used  in  con¬ 
junction  with  equation  (14)  to  obtain 


V  =  BM  (ns'  ) 
LN 


r  f 1 
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J0 F3(W) 

dw 

f  f2  (w, 

A  dw 
ns  J 
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Jo  x 

— 
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where  8  =  /XL 

M  (ns*  )  is  the  measured  "slope"  function  and  the  quantity 

in  brackets  is  a  known  function  that  can  be  calculated  numerically 
with  the  aid  of  a  computer.  It  has  been  calculated  for  the  case 


»’r  '  °-43 
D's  ’  4'85 

and 

L  =  1.6  Km 
For  this  case 

V  =  2.725  x  10~2  M  (0.43)  m  sec-1  (29) 

IN 

The  proportionality  constant  is  a  function  of  the  detector  spacing  and 
size,  the  source  size  and  the  path  length.  The  detector  spacing (s) 
and  size  are  known  and  the  path  length  will  be  known  even  in  the  field 
application,  since  a  range  finder  will  either  be  integrated  into  the 
system  or  coupled  in  externally.  However,  the  source  size  will  in 
general  be  a  function  of  turbulence  level  and  its  effect  on  accuracy 
and  calibration  needs  to  be  evaluated. 

The  instrument  generated  binary  data  tape  begins  with  the  follow¬ 
ing  sequence  of  symbols  which  indicate  the  beginning  of  a  data  block: 
1,0,//  of  bytes  in  data  record  (low  byte  first),  0,0.  This  is  followed 
by  the  data  block  which  consists  of  the  200  dita  words  from  the  array 
plus  7  extra  channels  that  can  be  used  to  record  in  situ  measurements 
of  wind  speed,  Cn2  etc.  for  evaluation.  A  data  word  consists  of  two  6 
bit  bytes  (12  binary  digits)  each  of  which  has  a  parity  (odd)  check. 
The  binary,  12  bit  analog  to  digital  conversion  gives  a  resolution  of 
one  part  in  4096  which  should  be  down  into  the  noise  level. 

The  data  is  sequenced  onto  the  tape  alternately  from  array  one 
and  array  two.  For  example,  1(1, 1,1)  is  punched  followed  by  1(1, 1,2), 
followed  by  I (1,2,1)  etc.  The  array  location  is  sequenced  along  each 
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row  to  the  end  and  then  shifts  down  to  the  next  lower  row  without 
retracing.  For  example  1(1,10,2)  is  followed  by  1(2,10,1). 

The  binary  data  on  the  instrument  generated  tape  is  converted 
using  the  Graduate  Center  PDP  11  computer  to  ASC  II,  and  a  new  tape 
generated  for  use  in  processing  the  data.  The  data  is  processed 
using  the  program  shown  in  Appendix  A  which  also  contains  a  sample  of 
the  computer  generated  output.  The  program  prints  out  the  data  from 
the  two  arrays  in  the  proper  sequence  and  then  calculates  and  prints 
out  the  average  intensity,  variance,  covariance  function  and  "slope" 
function.  The  data  used  in  this  example  is  hypothetical  and  is  used 
merely  to  check  out  the  program. 
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EXPERIMENTAL  PROGRAM 

Measurements  will  be  made  at  an  established  field  site  over  a  1.6 
Km  path  with  very  uniform  terrain.  A  small  building  complete  with  heat, 
light,  power,  communications  and  an  optical  table  completely  isolated 
from  the  building  was  constructed  to  house  the  experimental  system. 
Figure  7  shows  the  field  site  facility. 

Emphasis  during  preliminary  phases  of  the  experimental  work  will 
be  centered  around  measuring  the  covariance  function  and  average  cross- 
wind  velocity  with  the  experimental  system  and  comparing  the  results 
with  the  calculated  covariance  function  and  in  situ  wind  measurements. 

As  noted  earlier,  the  covariance  function  changes  with  source  size. 

In  order  to  control  this  parameter,  the  target  size  will  at  first  be 
restricted  to  a  diameter  of  4.85  /XL  and  the  beam  adjusted  at  very 
low  turbulence  levels  to  fill  the  entire  target.  The  calculated  co- 
variance  curve  for  this  case  is  shown  in  Figure  8  and  the  corresponding 
wind  weighting  function  is  shown  in  Figure  9.  It  can  be  seen  from 
Figure  10  that  the  weighting  is  peaked  up  toward  the  receiver. 

During  this  phase  of  the  work,  the  transmitter  was  operated  at 
the  field  site  and  photographs  were  made  of  the  received  scintillation 
patterns  from  the  double  pulse  ruby  laser.  Two  cameras,  an  electronic 
shutter,  a  rotating  disk  with  a  hole  and  mirror,  the  receiver  optics 
ana  electronics  to  synchronize  firing  of  the  laser  with  the  disk  were 
used.  The  rotating  disk  directed  the  scintillation  pattern  from  the 
first  laser  pulse  to  one  camera  and  the  second  pulse  to  the  other 
camera.  High  speed  infrared  film  (Kodak',  was  used  to  record  the 
patterns  with  the  Scotchlite  target  at  500  meters  and  1.6  Km.  Repre¬ 
sentative  scintillation  patterns  are  shown  in  Figures  10  -  12. 

Normal  looking  scintillation  patterns  were  obtained  at  both  ranges. 

The  primary  scale  size  appears  to  be  appropriate  for  the  path  length 
and  there  appears  to  be  no  discernible  speckle  present.  The  length 
of  the  longest  cross  hair  (from  cross  to  edge  of  aperture)  is  approxi¬ 
mately  8.5  inches.  The  patterns  were  recorded  at  night  under  very  low 
turbulence  conditions.  The  "rice  grain"  oattern  in  Figure  11  is 
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thought  to  be  due  Lo  a  very  light  rain.  As  can  be  seen  from  the 
photographs,  the  electronic  synchronization  with  the  disk  was  not 
perfect  and  consequently  part  of  the  patterns  were  occluded  by  the 
disk.  This  will  not  be  a  problem  with  the  actual  receiver  in  that 
a  "real  time"  response  will  be  available  while  adjusting  the  synchron¬ 
ization.  The  overall  results  from  the  photographs  of  the  scintillation 
patterns  tends  to  indicate  that  speckle  will  not  be  a  problem. 


29 


)0 


Figure  7.  Field  Site  Facility. 
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APPENDIX  A 


005  DIM  11(10,10),  12<10. 10),  L1<;10),  L2<10),  A<2).  V<2) 
006--D1M-X1.C  10,-10) ,  X2 < 10, 10  )  ,-Sl ( 10) >  Ml < 10 )  ,  Cl (2,  9), M3<9) 
010  OLD 

015R£M_R£AD  IN  CELL  DATA  AND  CONVERT  SEQUENCE 

020  FOR  J=© T09 

025  FOR  K=0TG9 

030  INPUT  I1(J, K), I2<J, K) 

0L2.5--NEX.T—K _ 

040  NEXT  J 

045  FOR  J=1TG9  STEP  2 
050  FOR  K=0TG9 
055  LET  M=9-K 
060  LET  L1<M)«I1< J, K) 

0-65— LE ,L_L2i.M ) .=  1.2_< J.,-K )  ...  _ 

0?@  NEXT  IC- 

dfS  FOR  K=0TO9 

080  LET  I1<J,K)=L1<K) 

085  LET  I2<J,  fO=L2<K) 

086  NEXT  K 

09J3_NEXT-_J_ _ _ _ 

095  REH_DATA  NOW  IN  NORMAL  MATRIX  FORM 
096  PRINT 

100  PRINT  "FIRST  ARRAY" 

101  FOR  Js0TO9 
110  FOR  K=0TO9 

115.  PRINT.  .m.J,.k),  ... 

120  NEXT  K 
125  NEXT  & 

130  PRINT  "SECOND  ARRAY" 

135  FOR  J=0T09 
140  FOR  f(=0T 09 
145  PRINT  I..2_<J ,  K), 

150  NEXT  K 

155  NEXT  J 

156  PRINT 

160  REM.CALCULATE  AVERAGE  INTENSITY  EACH  TUBE 
165  LET  S=0. 

170  LET  M=0. 

175  FOR  J=8 T09 
180  FOR  R=8T09 
185  L E T  S=S+I1‘J,K) 

198  LET  M=M+I2<.  J,  K> 

195  NEXT!; 

200  NEXT  J 

205  lET  A*.  1)=S/1@8 

210  LET  A; 2 ,=M/10C 
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215  REM  CALCULATE  VARIANCE 
220  LET  5=0. 

225  LET  M=0. 


250...FC<R-jL=0.T.Q5L _ — 

225  FORK=0TO9  --  ■  f 

240  LET  5=S+I1<J,  lO*ll<J>.K)  -  - - j - 

245  LET  M=M+I2< J> K)*I2< J>  fO  1 

250  NEXTK  - 

255  N£ja_ i  -  —  - - 

260  LET  V<1;=S/100.  -A<l>**<i> 

265  LET  V<2)=M/100.  -%< 2 >*&<%) 

270  PRINT  "PULSE",  “RV&Mli  INTENS1  TVS  “VARIANCE" 

2? 5  PRINT  1,  Ra),,  V<1‘> 

280  PRINT  2,  A<2)/>.Vt2) 

281  PRINT 

285  REM  CONVERT  1  DATA  TO  ZERO . MEAN  AND  NORMALIZE  SV  VARIANCE 

290  FOR  J=0TG9 

295  FOR  K=0TO9  ..  _ 

200  LET  Xl<J,k)»<U<J,KJ-A<l)>/SGR<V<l>> 

205  LET  X2<J,  K)-<12(J>  K)~R<2) )/SQR<V<2) ) 

210  NEXT  K 

315_.NE.X1-J _ 

220  REM  CALCULATE  COVARIANCE  FOR  9  SPftCINGS 
225  FOR  S=0TO8 
220  LET  M=10-S-2 
225  LET  52=0. 

240  LET  M2=0. 

245  FQR_J=0XQS _ _ 

250  LET  S1<J)=0. 

255  LET  Ml< J)=0. 


260 
265 
370 
37.5 
2  76 
280 
285 
290 
295 
400 
410 
415 
420 
425 


FORK=0TOM 

LET  Sl< J)=S1<J)+X1<J,  K)*Xi<J,  K+S+l) 

LET  Ml (  J  )  =M1  <  J  )  +X2( J >  KMX2CJ,  K+S+l) 

LET  S2^S2+51<jT*  . .  . ~ 

LET  M2=M2+M1< J) 

NEXT  J 

LET  Cl<li  S;=S2V10.  /< M+l .) 

LET  Cl<2i  S)=M2X10.  /<M+1) 

NEXT  5... _  ..  _  _ _ 

PRINT  "STEP", "COVARIANCE  PULSE  ONE", "COVARIANCE  PULSE 
FOR  J=0TO8 

PRINT  J+l,  C1<1  ,J),,  Cl <2,  J) 

NEXT  J 


TWO" 


426  PRINT 


420  REM  CALCULATE  SLOPE 
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425  FOR  S=0TOS 
4.40  LET  M=10.  -S-2 
445  LET  £2=0. 

450  FOR  J =0X09 
455  LET  Sl<  J>=8. 

A£euums£ma _ _  . .  . 

465  LET  S1<J)=S1<J)+X2< J,  K+S+1)*X1< J,  K)-X1<J* K+S+1)*X2< J* K> 

4?0  NEXT  K 

4?5  LET  S2=S2+S1<J> 

480  NEXT  J 

485  LET  M3<S>=1000.  *52/10.  /<M+1.  ) 

480  NEXT  S _  _  _  _ 

495  PRINT  "STEP'S  "SLOPE  OF  COVARIANCE" 

500  FOR  O'=0TO8 
510  PRINT  J+1,M2<J> 

520  NEXT  J 
520  STOP 

JWflLENfi. _ _ _ 


FIRST  ARRRV 

12  1 
2  12 


1  2  1 
2  12 


1 

2 

1 

2 

1 

2 

1 

1 

2 

1 

v*. 

1 

2 

1 

cL‘ 

i 

u. 

1 


1 

2 

1 

1 

‘I 

CL 

1 

1 

1 

U. 

1 

IM 

1 

1 

.V* 

1 
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SECOND  RRRftV 


H*-  K>  to  K  to  hi.  to  hi.  to  hi  to  hi  to  hi 


PULSE  AVERAGE  INTENSITY  VARIANCE 

. 1 .  . 1,5...  .25 

2  ^  1.  5  .  25 


STEP  COVARIANCE  PULSE  ONE .  COVARIANCE  PULSE  TWO 


.  „1 _  -1  *1 

2  1  1 

-.-2_  -1  -I 

4  1  1 

5  -1  -1 

6  1  1 

?  -1  -1 

8  1  1 

S  -1  -1 


STOP  AT  LINE  530 

READY _ 

READY 

SAVE 
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